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IN COMPRESSIBLE GAS STRFAM WITH SUBSONIC VELOCITIES* 

By L. G. Loitsianskii 


In the computation of the impellers of turbomachines, the Joukowski 
formula is applied; according to the formula, the lift force of the blade 
is equal to the product of the density of the gas, the mean vector of the 
velocities ahead of and behind the cascade, and the circulation about the 
blade . 

The presence of high oncoming relative velocities of the gas flow at 
the blade make it necessary to account for the effect of the compressi- 
bility of the air on the lift force of the blade. 

In 1934, Keldysh and Frankl (reference l) showed that, in the case 
of the isolated wing profile, the Joukowski formula retains its usual 
form for both incompressible and compressible flow. As far as is known 
to this author, no generalizations have yet been made to the case of the 
airfoil in cascade. ^ 

It is shown herein that, for subcritical values of the Mach number, 
the Joukowski formula may, with an approximation entirely sufficient for 
practical purposes, also be applied in its usual form to a cascade in a 
compressible flow, provided that the density is taken as the arithmetical 
mean of the densities of the gas at infinity ahead of and behind the 
cascade . 


This simple result may be useful for the computation of cascades, 
especially in the solution of the problem of the resistance of a cascade 
when it is necessary to separate the lift force from the total force deter- 
mined by the momentum theorem (reference 2). 


*"Obobshchenie Formuly Zhukovskogo na Sluchai Prof ilia v Reshetke 
Obtekaemoi Szhimaemym Gazom pri Dozvukovykh Skorostiakh . " Prikladnaia 
Matematika i Mekhanika. Vol. XIII, No. 2, 1949, pp. 209-216. 

■'"Very recently, there had come to this author's attention the as yet 
unpublished generalization of E. M. Berson, which differs from the one 
given herein. 
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1. Vector equation of lift force of an airfoil of a two-dimensional 
infinite cascade in an ideal incompressible fluid . - The flow about one 
of the airfoils of a cascade situated in a tube of flow having at infin- 
ity ahead of and behind the cascade the dimension equal to the pitch t 
in the direction of the axis of the cascade is considered (fig. l) . The 
density of the fluid is denoted by p, the velocities and pressures at 
infinity ahead of and behind the cascade by Wq, w 2 and pq, p 2 , 
respectively. When the momentum theorem is applied to the mass of fluid 
enclosed in the flow tube between the infinitely removed sections aq 
and eg , the equation for the vector lift force may be written as 


R = (Pq-P 2 ) t + p (t * Wp) Wp - p (t • w 2 ) w 2 (1.1) 


where the vector t is equal in magnitude to the pitch and is directed 
at right angles to the axis of the cascade downstream of the flow. 

The difference in pressures Pp~P 2 i- n equation (l.l) may be 
expressed in terms of the velocity vectors at infinity by the equation 

Pi - p 2 - f K 2 - w i 2) = I (V*i> * (V^i) d- 2 ) 

The mean vector velocity w m and the velocity w d characterizing 
the vector change of velocity of the flow in passing through the cascade 
are then introduced: 


= \ ('Wp+Wg ) 


w d = w 2 - -Wq 


Equation (1.2) then becomes 


(1.3) 


p - p 2 = pw m • w d 


(1.4) 


The condition of the conservation of mass along the tube of flow 


gives 
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(1.6), 

the 

vector 


( 1 - 5 ) 


( 1 . 6 ) 


plane and is directed perpendicular to w m toward the side determined 
by the vector product on the right side of equation (1.6). 


The magnitude of the vector R is given by 


R = pw m v d t 


(1.7) 


2 . Possible transformation of Bernoulli formula for compressible 
gas at subsonic velocities . - In the flow of an ideal compressible gas 
about a two-dimensional cascade, the condition of conservation of mass 
flow along a flow tube yields 


P 1 ^ ’ W l^ = P2 ^ (^ • l) 

and with the notation of equation (1.3), the equation for the vector 
lift force may be given in the form 

R=(P 1 -P 2 ) t+ Pl (t * w ] _) w ]L - p 2 (t • w 2 ) w 2 ■ (p x -p 2 ) t- Pl (t * w L ) w d 

( 2 . 2 ) 

where p^ and p 2 are the densities of the gas far ahead of and behind 
the cascade. 

The first component of equation (2.2) is now considered. When it 
is assumed that the motion is adiabatic, the known relations for the pres- 
sure and the density are 
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P = P 0 1 “ SI 


■ P ° 1 ' SI X2 + 2(^)2 - • • • ( 2 - 3 > 


k-1,2^ ‘ | n 1*2. 2-k *4 

p = p n l^ 1 - k^i x y =p ol 1 -¥s x + * 


2 (k+l) 2 


(2.4) 


where Pq and denote the pressure and the density, respectively, 

of gas adiabatically brought to rest, k is the ratio of specific heats 
of gas at constant pressure c and constant volume c^, and \ is the 

ratio of the velocity of flow w to the critical velocity of gas a + , 
expressed in terms of the velocity of sound a Q in the adiabatic stag- 
nation of the gas: 


\ = — 
a * 


[jL 

= V k+l 


2k £o 

a 0 = U k+1 Po 


(2.5) 


From expansion (2.3), the difference of the pressures becomes 

p i - p 2 - m p o (K 2 S -h 2 ) - l^Kiy (h 2+ ^2 3 ) + £ p 


(2.6) 


where 


p (\ 2 2 -\ L 2 )k/(k+l) 


/ k-1 

(l “ k+l 


(X 1 4 +W+*z 4 ) - • • • (2-7) 


and and are the values of the parameter X ahead of and behind 


the cascade. 
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From equation (2.4), the arithmetic mean p of the densities 
and pg can he expressed as 


p m = 2 (Pl + P2> “ P 0 1 “ 2(k+l) ^l^^ + 


where 


1 

C P " 2 


t - m + ("-si 


v 2 ' 
X 2 , 


,k-l 


- 1 + 


( 2 . 8 ) 


1 (\ 2 +\ 9 2 ) = (X U4) . . . . 


V 1 2 ' " 4(k + l) 2 V 1 


(2.9) 


When equations (2.6) and (2.8) are compared. 


k P 0 , 2 2,/ p m 

p i - ' EI ^5 <-2 -“i > ( ^ + £ P - £ p, 


=5 P W • V, 

K m m d 


i + TT (c.-0 

m 


(2.10) 


Use of equations (2.5) and (1.3) yield, respectively, 


k P 0 


1 ^0 


k+1 „ o 2 2 22 

p o a * P 0 a 0 


w_ 


2 1 - 
W 1 - 2 W m 


w. 


From equation (2.10) it follows that, with a certain error, the 
order of which will subsequently be discussed, the differences in pres- 
sures ahead of and behind the cascade may be determined by an equation 
similar to equation (1.4) for the incompressible fluid if the density of 
the gas is replaced by the arithmetic mean of the densities ahead of and 
behind the cascade. 
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3. Generalization of Joukowski formula to flow of a compressible 
gas at subsonic velocities about a cascade . - When equation (2.1) is 
used, the second component of equation (2.2) "becomes 


P i ( t * w l) w d = p m (t ' + p l (t ‘ W l } " p m^ ' V 


V- 


= P m ( t • w m) w d + 


P 1 P 2 _ 


p. ( t ' w m) " P m ( t • w m)| -d 


m 


rp l p 2 


i) 


= P m ( t * w m) w d + l— ' P m ) ( t • w J w d 


, Pm "PlP2 . 

= U - ~ 2 J Pj t ‘ w m) w d 


m 


( ^£^1 
y± + % 


P (t * 

I m m' d 


In order to estimate the second component in the brackets, note 
that, according to equation (2.4), 


P 1 ' p 2 - k+1 


P 1 + P 2 = 2p 0 


0 ^ 2 - X l 2) ” 2(T/t) (\ 2+X 2 2) + * * *] 

" 2 (k+1) ^1 +^ 2 ) + • • J 


(3.1) 


Hence, 


' P l~ P 2 | 

Pl+Pg/ 


1 -2 


- \ 2 ) 2 + (^l 2+ ^2 2 ) + - - -J (3-2) 


4(k+l) 

Therefore, 

Pl (t • w 1 )w d = p m (t • w m )w d - e 'p m (t • v m )w d 


(3.3) 
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where, from what has just been shown and when, for convenience, 
K = l/(k-l), 

fP 1 f p l~ p 2^ 


P, + 


p 2 ) 


~ 4 


m 


|l - (k-l)/(k+l)X 1 2 J K - ji - (k-l)/(k+l)\ 2 2 J 


12 


1 - (k-l)/(k+l)^ 1 2 * - 1 - (k-l)/(k+l)^ 2 * 


4(k+l) 


(x 2 -\ 2 ) 2 


2 v/v 2 ^1 


k-1 

1 + k+1 


(\i 2 +\ 2 2 ) 


(3.4) 


The equation of the order of the error e p ' will be ignored for the 

present; equation (2.2), on the basis of equations (2.10) and (3.3), may 
then be represented in the form 


R “ P m ( w m ’ w d> 1 - p m ( t * v m> w d + e R = p m w m X (* x ^ d ) + e R ( 3 - 5 ) 


where represents a certain small vector 


P0 


£r = T ( e p“ € p ) p m K* w d) t + £ p' p m (t * w m } w d 
m 


(3.6) 


the order of magnitude of which will subsequently be discussed. 

When this small vector is neglected, generalization of the 

Joukowski formula to the flow of a compressible gas about an airfoil in 
cascade yields 
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When equation (3.7) is compared with the previously derived equa- 
tion (1.6) for the flow of an incompressible fluid about a cascade, the 
following result is obtained: At subsonic velocities, the lift force of 

an airfoil in cascade can be approximately determined by the Joukowski 
formula for the incompressible fluid if the density of this fluid is 
equated to the arithmetic mean of the densities of the gas at a great 
distance ahead of and behind the cascade - 

The absolute magnitude and the direction of the. vector R with the 
assumed approximation will now be determined and the difference from the 
magnitude and the direction of the vector R determined for the incom- 
pressible fluid by equation (1.6) will be shown. 

The vector product t X w d is perpendicular to the plane of flow of 
the gas and therefore to the vector w , which yields 

R=p m w m| Ix ^d| ( 3 - 8 ) 


In contrast to the incompressible case, in the compressible case the 
vector w d is not perpendicular to the vector t, but is inclined by a 
certain small angle to the axis of the cascade (fig. 2). Hence, in con- 
trast to the condition of perpendicularity (equation (1.5)) in the case 
of the incompressible gas, equation (2.1) gives 


t • Wd = t • (wg-w d ) = t 


h p 2 n ' £ »in) * • "in) (t - Tj 


t . v m - 1 . (» 1 + v 2 ) -i-(~ pin + ^ p 2 n) - 1 (* • pj.V + \ 


Hence, 


2‘^d 2 

t • w m P d + Pg 


P 1 p 2 P 1 p 2 


m 


- z f7 


(3.9) 


Whence, from equation (3.2), 


t • w d = 2 


1 d,V)(t-;)|i. 


w m = k + 1 ^' v 2 ~' v l 


1 + • • • 


(3.10) 
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Denoting the angles between the vector t and the vectors v m and 
w d by © m and 9^ yields the following expression for the modulus 
t X w d |: 


t x w d | = tw d sin 9 d = ^|t 2 w d 2 - t 2 w d 2 cos 2 9 d » V t2 ^ ' ^ * w d )2 


- 


2 4 € / (t • wj 2 = tw d \j 1 - 4c, 




cos 2 9. 


m 


= tw- 


1 - 2c 1 (— ) cos 2 9 

P UJ m 


(3.11) 


The modulus of the vector R is thus given by 


R = Pm w m w d t 


1 - 2e n ' 


' W mY 2 


p 


cos 9, 


m 


p ¥ ¥.t 
K m m d 


(3.12) 


The direction of the vector R for compressible gas also remains 
perpendicular to the vector of the mean velocity w m (fig. 2). 

4. Estimate of accuracy of generalized Joukowski formula . - In order 
to answer the question of the possibility of application of the previously 
derived equations, it is necessary to estimate the order of the neglected 
terms and also the interval of Mach numbers in which, with a preassigned 
accuracy, the proposed approximate formula may be applied. 

For this purpose, consider the estimate of the order of magnitude of 
the vector e R . The simplest method of estimating is applied in writing 
the evident inequality 


P ° I I 

€R ~ ^ | e p " £ p| cos 9 md + VVVV cos 9 m 


or, when the cosines are replaced by their maximum values. 


£ R < 


fo 

Pm 


£ P " C P 


+ € ' ) p w w t 
P / m m d 


(4.1) 
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The ratio p /p is readily estimated from the known equations: 
0 m 


P 0 = p l 1 + ~ M ' 


k-1 2 

w ^ 

2 


2 


1 ) “ p 2 V 1 + 2 M 2 


K = 


k-1 


Setting M-j_ = M 2 = 1 yields 


p 0 < p 






k-1 


p o < p 2 \ 1 + ~r - p : 


c 


\ K 
k+1 \ 


2 V 2 / 


Assuming, however, that M-^ = Mg = 0 yields 


p 0 > P 1 P 0 > p 2 


The result immediately follows: 




i< ~ < 


m 


2 / 


K = 


k-1 


(4.2) 


The estimate of the remaining magnitudes entering the parentheses of 
equation (4.1) are now considered. For this purpose, the following 
simple device is used: The values k s= 4/3 = 1.33 and k = 3/2 = 1.5, 

between which the region of most practical values of k is found (for 
air and certain other gases), are considered and for these values of k. 


the values of €, 


and c. 


are computed. Similar computations may 


’P' P* ~P 

easily be made if desired for the interval 5/4 ^ k < 4/3, and so forth. 


For the chosen values of k, the exponents of the terms in equa- 
tions (2.7), (2.9), and (3.4) become integers and the entire computation 
is easily carried out. Then, for k = 4/3, 
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,, 1 X 2 n 4 1 X 2 \ 4 

(1 - 7 *1 ) ~ (1 - 7 *2 } 

y ( x 2 2 - ^i 2 ) 


1 + n <^i 2 + S 2 > 


s= ~ + ^ 2 ^2 2 + ^2^) —5 ( ^2 2 + ^2^ + ^2^ 


C P 2 


£-^i 2 ) 3 + ( 1 -7V) B ]- 1 + H<\ a+ S* ) 


“^2 (^l 4 + ^2 4 ) " (^i 6 +^ 2 6 ) 

2-7^ 2-7° 1 ^ 


1 ^ p l _p 2^ p 0 2 [/ 1 


£ '_ — - — _ 

P " 4 2 .2 

Pm 4 P m 


4^ 2 Uj ( ^ 2 ' X;:2) 


i- i^zY - h-i\A 


( x - 7 V*) (* - 7 ^2 2 ) + (* - 7 vj 
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On the basis of the obtained formulas, an estimate of the expression 



C R 

Pm w m w d t 


(4.5) 
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for any values of k In the interval 4/3 < k < 3/2 is given: From 
equations (5.2), (5.3) (Ed. note: equations (4.2) and (4.3)"1, and so 
forth, for k = 4 / 3 , 



or 

p 2 

0.033 (X-i 2 - \ 2 2 ) 2 <— l e p " c p| + € p ' < 0.133 (\ 2 - \ 2 2 ) (4.6) 

I I " 

and for k = 3/2, 



or 

0.032 (>^ 2 - X 2 2 ) 2 < ^ | e p "£ p | + e p ' < 0.108 (^2 - \ g 2) 2 ( 4 . 7 ) 

Thus, for the magnitude £„, the estimate is 

0.033 ( V ” X 2 2 ) 2 C p C / w t < 0 - 133 - X 2 2 ) 2 ( 4 . 8 ) 

m m d 

From the parameters and A- 2 , the Mach numbers M p = w p /a p an< ^ 

M 2 = Wg/ a.^ are obtained, where a p and a 2 are the velocities ahead of 
and behind the cascade. By the known equations of gas dynamics. 
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X 2 = fcfl M l 2 

1 2 1 + | (k-1) M x 2 

X 2 k+1 M 2 2 

2 2 1 + -| (k-1) M 2 2 

Whence, 

m 2 M 2 

X 2 - X 2 iti **1 ~ M 2 

1 2 [i + | (*-U «x 2 ] [i + -5 ( k -D «?j 

Hence, 

TZi (»1 2 - M 2 2 ) 2 < (kl 2 - k 2 2 ) 2 <( £ g i ) (M X 2 - M 2 2 ) 2 (4.9) 


Thus, in the interval 4/3 < k < 3/2, the following estimate may he 
used in place of equation (4.8): 

2 

0.03 (M-j 2 - M 2 2 ) 2 <— £_ <0.2 (M x 2 - Mg 2 ) (4.10) 

p m m m 


Equation (4.10) is a very rough estimate; in fact, the modulus of 
the vector e R will he much less. From this simple estimate, however, 

the possibility of practical application of the approximate equation (3.7) 
is evident. Thus, for example, if at the entrance to a turbine cascade 
the Mach number is equal to M-^ = 0.2 and at the exit Mg = 0.7, the 

relative error will not exceed 4 percent, even for the very rough estimate 
assumed. 


Translated by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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